#include "isotropic_atmosphere.h"
#include <math.h>

IsotropicAtmosphere::IsotropicAtmosphere(QString layerName,
					 double hBottom,
					 double hTop,
					 QObject *parent) :
    AtmosLayer(layerName, hBottom, hTop, parent)
{
}


void IsotropicAtmosphere::solarScatCoeff(int thetaS, int phiS,
		    linearPolar *lightOut)
{
    double lightIncident[3];
    lightIncident[0] = m_solarIntensity ;
    lightIncident[1] = 0;
    lightIncident[2] = 0;

    lightOut->i = lightIncident[I]/4;
    lightOut->q = 0;
    lightOut->u = 0;
}

void IsotropicAtmosphere::scatCoeff(int iThetaS, int iPhiS,
	       int iThetaI, int iPhiI,
	       double lightIncident[3],
	       linearPolar *lightOut)
{
    lightOut->i = lightIncident[I]/(4*M_PI);
    lightOut->q = 0;
    lightOut->u = 0;
}
